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Abstract 

We present a new recursion and Hamiltonian operators for the Viallet equation. This new recursion 
operator and the recursion operator found in fT| satisfy the elliptic curve equation associated with the 
Viallet equation. 

1 Introduction 

It is well known that integrable equations possess infinitely many generalised symmetries. These symme- 
tries can be generated by so-called recursion operators [2 [3], which map a symmetry to a new symmetry. 
For example, the famous Korteweg-de Vries (KdV) equation 

Ut = Uxxx + QuUx 

possesses a recursion operator 

^ = dI + 4u + 2uxD-^ , 

where stands for the right inverse of total derivative D^. This is a second order pseudo-differential 
operator which action is defined on elements of the space Im Dx ■ The infinite hierarchy of symmetries of 
the KdV equation is generated by 

ut^=W{ux), j = 0,l,2,---. 

Majority of the known recursion operators [3] are weakly nonlocal [5j (Some exceptional ones can be 
found in [S] and were further studied in [7|), that is, their nonlocal terms are a finite sum of terms of the 
form KD~^5uP, where K is a generalised symmetry and 6uP is the variational derivative of a conserved 
density p of the equation. For the KdV equation, any weakly nonlocal recursion operator is a polynomial 
of 3? with constant coefficients. However, this is not true for elliptic models such as the Krichever-Novikov 
equation and the Landau-Lifshitz equation. 

In [8], Demskoi and Sokolov showed that the commutative ring of all weakly nonlocal recursion 
operators for elliptic models is isomorphic to the coordinate ring of the elliptic curve. For example, the 
Krichever-Novikov equation 

3Uxx 

— ^xxx 7^ I ) 

^ Ux Ux 

where P(ii) is an arbitrary quartic polynomial of dependent variable u, possesses one weakly nonlocal 
recursion operator 3?iof order 4 and another one of order 6. These two operators are related by the 
algebraic equation corresponding the elliptic curve associated with the Krichever-Novikov equation: 

^1 = ^1- - e, (1) 
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where (p and 6 are proportional to two modular invariants of a quartic polynomial P(u). 

In this letter, we first present a new recursion operator valid for both the Viallet equation and 
Yamilov's discretisation of the Krichever-Novikov equation. We state that the new recursion operator 
and the operator proposed in [UIU] satisfy to the algebraic equation similar to ([T]). 



2 The Viallet equation 

In this section, we give necessary definitions and a short account of relevant results on our recent study 
of the symmetries, conservation laws and recursion operator for the Viallet equation. 
The Viallet equation is a difference equation on 1? given by 

Q := aiUo,oMi,oWo,iWi,i 

+02 (■"0,0^^1,0^^0,1 + ^^1,0^^0,1^^1,1 + Uo^lUi^iUofi + Ui^iUofiUifi) 

+a3(^io,o^^l,o + ^^0,1^^1,1) + a4(ui,oUo,i + ^^0,0^1,1) 

+05 (■"0,0^^0,1 + ^^1,0^^1,1) + a6('"-o,o + uifi + no,i + ni,i) + 07 = , (2) 

where Oj are arbitrary complex parameters such that the polynomial Q is irreducible. Here the dependent 
variable u is a complex- valued function of independent variables n,m £ We have two commuting shift 
maps S and T defined as 

5 : M i-> ni^o = u{n + 1, m), 
T : u ^ uq^i = u{n, m + 1), 
S^T'^ : u I-)- Up^q = u{n + p,m + q). 

For uniformity of the notation, it is convenient to denote the "unshifted" function u as iio,o- 

Equation ([2]) was identified by Viallet with the vanishing of its algebraic entropy [10]. By a point 
fractional-linear transformation, it can be reduced to Adler's equation, also referred as the Q4 equation 
in the ABS classification jllj . In fact, all of the ABS equations can be obtained from the Viallet equation 
by a simple specialisation of parameters. 

With affine linear equation ([2]), following [12], we associate a bi-quadratic polynomial h (the discrim- 
inant of Q) defined by 

h{uofi,uifi) = Q<9no,i^«i,iQ ~ duo^iQdu^^Q (3) 
and a quartic polynomial / of no,o given by 

AfRne-linear equation ([2|) is covariant with respect to the Mobius transformations 

aUn,ra + f3 . a 

Classical modular invariants §2,93 £ C are defined as (see |12] , [13] ) : 

92 = (2//^^ - 2/7'" + (/")') , (4) 
53 = ^ (l2//'7'^ - 9{f'fr - (if iff + Gf'ff'" - 2{f"f) . (5) 
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Thus with a difference equation defined by an affine-Unear polynomial Q one can associate a plane 
algebraic curve 

in the Weierstrass form. If the modular discriminant of the curve is non-zero §2 — 27g1 then the 
curve is elliptic \T3\. 

Following [12j we introduce the relative invariants 12,^3 G C: 

h = l {hdldl ^h - {duh){dudl^^^h) - idu,,Mdldu,^oh) + {dlh){dl^,h)) (7) 
+^{dudu,^,h)\ (8) 

I ( ^ dlh \ 

/3 = -det du^^h dujoduh du^^dlh , (9) 

We note that the relative and modular invariants are not independent. 
Proposition 1. There is a syzygy among 12,13,92 o,nd g^, namely 

II = All - 92I2 - 93 ■ (10) 

Proof: The statement of this proposition can be verified by a direct computation. ■ 
The Viallet equation ([2]) is invariant under the involution 

ui,o ^ uo,i, as ^ 05. (11) 

This property enables us to study symmetries, recursion operators and conservation laws only for the 
one direction of the lattice and to recover the complementary set for the other direction. We denote J^s 
the field of rational functions of variables {unfi \ n € Z}. It is a difference field with automorphism S. 

Definition 1. Assume K ^ Fs depends on a finite set of variables {un,o \ ni < n < 712} with du„_^ qK 7^ 
and du„,^ „K ^ 0. We say K is a symmetry of equation (0) of order (ni,n2) if 

Dq{K) = for all solutions of Q = 0. 

Here Dq is the Frechet derivative of Q defined as 

n, = ^ST^^T+^S^^. (12) 

Equation ([2]) possesses a generalised symmetry of order (—1, 1) [HI [T5l [T]: 

^(1) = h h_ (13) 

The corresponding generalised symmetry for equation ([2]) in T direction is 



h 1 



uo,i - Uo-l 



-B h 
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with the polynomial h = Q du^ ^^du^^Q — du^^Q du^^Q . It is related with (fT3]l by the involution (fTT]l . 

For convenience from now on, we will use one-index notation to denote shifts of elements J-g by S, 
e.g. K^^^ = and K = 5"(/i(no,o, tii,o))- We shah omit the index zero for unshifted functions. 

Often we shall use notation 



w = and Wk = S w. (14) 

In this notations a generalised symmetries of the Viallet equation of order (—2, 2) can be written as 

[a da ESI [I]: 

= h h_iw^{wi + w-i) . (15) 

In our recent work [1, ,9], we proved that infinitely many generalised symmetries K^"'\ ovdK^"'^ = 
(— n, n) of the Viallet equation ([2]) can be generated by a recursion operator 9\ 

J^{2n-l) ^ 9^n-l(i^(l))^ |^(2n) ^ (if (2) ) ^ ^eN. 

Here the recursion operator is a second order pseudo-difference operator of the form 

Is' + h h^, wW^S'^ + 2^(1)^(2) ^ 1^ 

—up' {h-i hiw\ + h-2 huPiP) + 

+ 2 (5 - l)-i o f - ^ ) + 2 (5 _ o I ^ _ fi^ I . (16) 



9^ = hh^iuPw 





For instance, lH(/ir*^^)) gives rise of the symmetry of order (—3,3) for equation ^ [Ij as follows: 

if (3) ^ ^(^(D) ^ ;^^_^^2 J^^(l)^2 ^ if ^ ^ if (Dif (2) . (17) 

Operator 9^ is a sum of a difference operator with coefficients in Ts and a non-local pseudo-difference 
operator. The nonlocal part of the operator ()16p is a finite sum of the form a (5 — 1)~^ o 6 with 
a, 6 G J-'g. Such operators belong to the type of weakly nonlocal pseudo-difference operators [Il[9]. It is a 
difference analogue of weakly nonlocal pseudo-differential operators studied by Maltsev and Novikov [5]. 

Recursion operator 9^ can be presented as the product y{ = T-Lol of the Hamiltonian operator T-L and 
the symplectic operator I, where 

n = h.ihhiw'^wlS - S-^h^ihhiw''wl+ 2K^^^S -ly^SoK^'^ (18) 
+ 2i^(2) (5_i)-ioA'«, 

T = \s-S-^\. (19) 

The Hamiltonian operator % maps the variational derivatives of conserved densities to symmetries while 
the symplectic operator X maps symmetries to the variational derivatives of conserved densities [lll^. 



3 A new recursion operator for the Viallet equation 

In this section, we present new weakly nonlocal recursion and Hamiltonian operators for the Viallet 
equation ([2]). The new recursion operator and the recursion operator ()16p satisfy to the agebraic equation 
corresponding to the elliptic curve associated with the Viallet equation. It is a product of the new 
Hamiltonian operator and the symplectic operator I given by (jlOp . 
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Theorem 1. A weakly nonlocal pseudo- difference operator 



m = —S^ H 5^ + I — h 

h2 hi \ h h 



hi I V 



+ 



/l \ /l_2 / h^2 h-3 



+ — 


h-i 






h-i 


h 






h-i 


h 



5-1 + ^5-2 + 



-3 



+2 /fW (5 - 1)-' o I 11^ - iii- 1 , (20) 

where Wi is defined by ^141); h is the discriminant of Q given by K^^\j = 1,2,3 are symetries of 
equation (0j given by i TOj) . I[T5\) and (TT^ , and 

A = h-ihhih2up'wlwl ; (21) 
B = 2h^iw^l^hiwl{K^^^W-i-]^wdu,^,h + -^du,,Ai,oh) + K^^^ , (22) 

is a recursion operator for equation 

In the proof that y{ p6p is a recursion operator of equation cf. Theorem 1 in [9], the key 
observation was that there exists a constant fi 

= (a4 — a3)(a3aia7 — 0305 — 20702 + 4050205 ~ 2oiOq + 04O1O7 — 04O5) , 

such that the operator d\ — fj, can be represented in the factorised form 

= Mf^S+ 



duifi duofi 

where is a weakly non-local pseudo-difference operator, explicitly given in Lemma 2 in In the case 
of the operator ()20p one can check that for 

jl = h + ^/i(oi07 + a\ + a\ + a\- 4o206) 

the operator tiK — fi also can be represented in the factorised form 

\auifl duofl 

with a certain weakly non-local pseudo-difference operator M (the explicit form of M is rather long and 
we omit it in this short paper). The rest of the proof of this theorem is not different from the proof of 
Theorem 1 in [S] which states that 9^ is a recursion operator for the Viallet equation. 

Proposition 2. The recursion operator D\ can be written as the product = % oX of the Hamiltonian 
operator % and the symplectic operator X given by iiy\} . where 

n = AS^ - S-^ A + BS-S-^B + K'^^^iS-l)-^ {S + 1)0 k'^^^ 

+ 2 K^^^ (5 - l)-^S o + 2 (5 _ o ^-(1) , (23) 

and A and B are defined in [2l\) and 
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Computing the product T-LoX one can directly verify that the result coincides with ()20p . Operator 
% is obviously skew-symmetric. The proof of the Jacobi identity for the Poisson bracket corresponding 
to T-i is similar to one for the operator % (Proposition 2 in [9]). 

The central result of this paper is formulated in the following Theorem. 

Theorem 2. (i) The recursion operators 9^ / lid]) and !i2U\) satisfy the algebraic equation 

(29^ - hf = 4(9^ + hf - 52(5^ + h) - 93, (24) 

where h,h,g2 oind are the relative and modular invariants (0), (0), 0j and ([^. 
(ii) Operators 9i and 9^ commute. 

Every pseudo-difference operator can be uniquely represented by its formal Laurent series and there 
is a unique formal series corresponding to a formal inverse of the operator [1]. For example operator 9^ 
can be represented by 



and its formal inverse by 



5^-1 ^ I ^-2 '2.K^_Ik[_1 ^ 

It follows from Theorem 1, that equation ([2]) has the first order formal recursion operator 

m = d\L°^l^ = h^iw^S + • • • . 

Logarithmic residue />o = res log = log(/i_itf^) and = res9^'^', k = 1,2,3,... are densities of 
canonical conservation laws for the difference equation ([2]) (see Theorem 4 in [Ij). 

In [J we claimed that we can compute a square root of the recursion operator 9^ ()16p . although for 
a general second order formal series a square root may not exist. 

Proposition 3. There exist first order formal series rQ + riS ^ + --- and R = h^iuP'S + 

fo + riS^^ + ■ ■ ■ with coefficients rk,fk € Ts such that D\l = R^ and 9^i = R^. 



Proof: The existence of the square root of 9^ ()16p and a cubic root of 9^ (I20p follow from Theorem 
[2] and the following obvious Lemma: 

Lemma 1. Let A and B be two commuting formal series and ord^'^ > ord^B for an integrer k > 0. 
Then the k-th root C = {A'' + B)^ is given by the formal series C = A + J2j>o^i^^ j)^'' •^^''''^ > where 
C(^; j) = j! 0^=1(1 — i + 1) are generalised binomial coefficients. 
Indeed, after re-arranging equation ([M]) as 

9^ = ^(2^ - /3)2(9^ + hr^ -h + j(^ + h)-^ + j(9^ + h)-^ 
we can apply Lemma [T] with 

k = 2, A = ]^{2^L-hmL+hr\ B = -h + ^{^L+hr^ + ^{^L+hr^. 

We have ord^ = 1, ord^B < and [91,9^] = implies [A,B] = 0. Then i? = C is a first order formal 
recursion operator with coefficients in J^s, such that 9\l = R^- In a similar way one can find R, such 
that D\l = R^- ■ 
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4 Yamilov's discretisation of the Krichever-Novikov equation 



As pointed in [9j, the operator (|16|) is a recursion operator of differential-difference equation ut^ = K^^\ 
wfiicli can be identified as Yamilov's discretization of the Krichever-Novikov equation (YdKN) [17], cf. 
equation (V4) when = in [TH]. Indeed 

u,, = = (25) 
«i,o - U-1,0 

where 

R{uifi, u, ii_i,o) = {hu^ + 2b2U + 63)iii_o'U-i,o + {hu^ + hu + ^5)(^i,o + ^^-1,0) + hu'^ + 265"" + h . 
and 

bi = aios - a|; 62 = |(a2a3 + ai«6 — 0205 — 0204); ^3 = 0206 — 0405; 

&4 = ^(ai + ai«7 - a| - Og); 65 = ^(0306-^0207-0406-0506); 66 = 0307-05. 
It is straightforward to check that 

h{u, uifi) = R{uifi, u, uifi). (26) 

Using (p5]) . ([26|) we can express *H (fT6|) . H (fT8]) and I (fT9|) in terms of the polynomial i?(ni^O; ^-i,o)- 
In [9] it has been shown that H, I and 9\ are Hamiltonian, symplectic and recursion operators for the 
YdKN equation respectively. By direct computation, we can show that 

Lj^(i)'H = 0, 

where is the Lie derivative in the direction of the vector field K^^^ (for a detail definition of the Lie 

derivative see [9J). Thus ?^ is a Hamiltonian operator for equation ()25p and 9^ is a recursion operator of 
YdKN (|25p . The mentioned above canonical conserved densities pk, k = 0,1,2, .. . are densities of local 
conservation laws of the YdKN equation. 
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